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Abstract. The Weibel instability is investigated with PIC simulations of an initially 
unmagnctized and spatially uniform electron plasma. This instability, which is driven 
by the thermally anisotropic electron distribution, generates electromagnetic waves 
with wave vectors perpendicular to the direction of the higher temperature. Two 
simulations are performed: A 2D simulation, with a simulation plane that includes 
the direction of higher temperature, demonstrates that the wave spectrum is initially 
confined to one dimension. The electric field components in the simulation plane 
generated by the instability equalize at the end of the simulation through a secondary 
instability. A ID PIC simulation with a high resolution, where the simulation box 
is aligned with the wave vectors of the growing waves, reveals details of the electron 
phase space distribution and permits a comparison of the magnetic and electric fields 
when the instability saturates. It is shown that the electrostatic field is driven by the 
magnetic pressure gradient and that it and the magnetic field redistribute the electrons 
in space. 
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1. Introduction 

The filamentation instability or beam-Weibel instability has been proposed as a 
mechanism to magnetize the early universe [1-6] , to provide the strong magnetic fields 
for the afterglow emission of gamma-ray bursts [7-10] and supernova explosions [11] and 
to explain the heating processes in the pulsar winds [12, 13]. It is also important for 
the inertial confinement fusion in laser-plasmas [14]. The filamentation instability has 
been investigated analytically [7,15-17], numerically [18,19] and also in a laboratory 
experiment [11], where it has been discussed as the precursor to the formation of 
astrophysical shocks of gamma-ray bursts and supernovae. 

However, the filamentation instability is not the only instability that can generate 
magnetic fields. Here, we focus on a second one, the thermal anisotropy-driven Weibel 
instability (TAWI), which we study with particle-in-cell simulations. There are several 
analytical and numerical works of the TAWI, e. g. [13, 20-24]. We illustrate the physical 
mechanism behind the TAWI for a plasma, in which the temperature along one direction 
exceeds significantly those along the other two directions. The typical microcurrents of 
the electrons are highest along the direction with the high temperature. The magnetic 
interaction force that causes the deflection of electrons is thus strong orthogonal to this 
direction. At the same time, the typical particle's thermal energy orthogonal to the 
direction with the high temperature is relatively weak. The particle's thermal motion 
can, in this case, not oppose the structure formation. Current channels form with an 
axis, which is aligned with the direction of the high plasma temperature. 

The original work of Weibel [25] considers a bi-Maxwellian distribution with one 
component of lower temperature yielding a one dimensional wave growth with wave 
vectors k in the direction of lower temperature. In contrast, we choose only one hot 
direction. The wave growth will then occur in the two dimensional plane perpendicular 
to that direction and give rise to the generation of a magnetic fleld in this plane. 

The temperature anisotropy we consider here can be the result of a Buneman 
instability [26], which emerges for example at the terrestrial foreshock [27]. The latter 
is a result of the backstreaming of solar wind particles, which do not pass the terrestrial 
bow shock and move back into the solar wind domain. The backstreaming particles 
form a foreshock upstream of the shock. A Buneman instability develops between the 
shock-reflected ions and the electrons of the solar wind, if we assume for simplicity that 
no magnetic fleld is present. The phase velocity of the waves is comparable to the ion 
beam velocity [28], which we denote here as vj, in the electron rest frame. The electric 
fleld amplitude grows and the electrons are trapped in phase space holes, when the 
kinetic energy of the electrons in the wave frame is not sufficient to escape the wave 
potential [29, 30]. They are accelerated along the beam velocity vector to a few times the 
beam speed. Since the spectrum of the waves driven by the Buneman instability is not 
strictly mono-directional, the electrons also experience a, usually weaker, acceleration 
away from the beam velocity vector. After the heating of the electrons, the Buneman 
instability is quenched. The electrons are heated more in the ion beam direction than 
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in the plane orthogonal to it and a temperature anisotropy has developed. 

The proximity of the terrestrial foreshock permits a direct access by in-situ 
observations. Therefore the particle configurations are well-known and the ion beam 
velocity could be determined to some hundred km/s [27], which is comparable to the 
speed of the Solar wind in the reference frame of the bow shock. We assume that this 
physical scenario can be applied to the foreshocks of supernova remnants (SNR), where 
the bulk flow speed of the supernova blast shell and thus of the SNR shock has been 
determined to be less or equal to 0.2c [31]. The ion beam speed should be similar to 
the velocity of the shock in the upstream reference frame. 

Realistic values for the interstellar medium, into which a SNR shock is expanding, 
would be BjsM = 1 nT and uism ~ 1 cm~^ [33] , which would give a ratio of the plasma 
frequency to the electron cyclotron frequency of ~ 300. Thus, we can neglect the 
magnetic fields. Experimental evidence indicates that the magnetic field close to SNR 
shocks is amplified to values ~ 100 nT, presumably by the shock-accelerated cosmic 
ray particles [34,35]. Recent PIC simulations [36] seem to suggest though, that the 
magnetic field cannot be amplified by the cosmic rays to values much higher than the 
amplitude of the background magnetic field. 

The TAWI we consider here is driven by the electron heating by the shock-reflected 
ion beam, which we find even at unmagnctized shocks [37]. It can thus amplify a 
magnetic field from fiuctuation levels and constitute the first stage of a chain of magnetic 
field amplification mechanisms upstream of a SNR shock. Our PIC simulations evidence 
that the peak magnetic field, which results from the TAWI, would be 8 nT for a plasma 
density of 1 cm~^, well in excess of the ISM magnetic field. These 8 nT constitute a 
minimum estimate, because the electron density, which scales this peak field, is higher 
in the foreshock than in the ISM. 

We perform two simulations and compare them to [22]: First we choose a 2D 
simulation plane that includes the direction of higher temperature, which we refer to as 
the parallel direction, and one perpendicular direction with a lower temperature. We 
show that the wave spectrum is one-dimensional during the early stage of the instability. 
In reality, this wave spectrum would be spread out over the plane that is orthogonal to 
the parallel direction. The dimensional reduction is, however, beneficial for the study of 
the interplay of the electromagnetic fields and the key plasma structures. We discuss the 
temporal behaviour of the electric and magnetic energy densities and analyse the field 
data, demonstrating that only that component of the magnetic field is amplified, which 
points out of the simulation plane. The magnetic pressure gradient generates initially 
an electric field in the perpendicular direction. A secondary instability transfers electric 
field energy from the perpendicular into the parallel direction. 

Second we perform a ID simulation along the perpendicular direction. We 
demonstrate that a weak electric field in the parallel direction is generated by Ampere's 
law. The connection between the much stronger electric field in the perpendicular 
direction and the magnetic pressure gradient is demonstrated also in this setting. The 
growing magnetic field pushes the electrons together, yielding an increasingly complex 
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layered structure in the velocity distribution. 

The paper is structured as follows. The initial conditions and the simulation setup 
are given in the section 2. The results of the 2D and ID simulations are presented in 
the sections 3 and 4, respectively. Section 5 discusses the results and the future work. 

2. The instability, the initial conditions and the simulation method 
2.1. The linear instability 

We discuss the TAWI for electrons in a homogeneous, coUisionless plasma and leave 
the ions as neutralising immobile background, which is equivalent to an infinitely large 
ion mass (mj — > oo). The plasma is initially unmagnetized {Bq — 0). A temperature 
anisotropy A — {vth\\/vth±y — 1 > in the electron particle distribution 



leads to growing electromagnetic oscillations with a wave vector k perpendicular to the 
direction of higher temperature. The velocities are normahzed to the speed of hght c. 

We denote the direction of high temperature as the x-direction {vth\\ = Vthx)- Hence, 
the wave growth of the TAWI mode can occur in the y-z plane. We can limit our analysis 
to the y-dircction due to the isotropy in this plane {vthL = Vthy), accordingly Vth± < Vth\\. 

The chosen thermal velocities arc Vth± = 8.83 ■ 10^^ and Vth\\ = 8.83 ■ 10^^, yielding 
A = 99. The parallel thermal velocity is a few percent of the speed of light, which 
matches with the ion beam speed upstream of the SNR shocks. The perpendicular 
electron thermal velocity corresponds to a thermal energy of 40 eV. 

The dispersion relation for the linear phase of the instability [22, 38] 



describes the growth of the electromagnetic oscillations with a wave number k in units 
of the inverse electron skin depth ujp/c and its associated linear growth rate a in units of 
LJp. The frequency is purely imaginary (a; = za), thus the oscillations do not propagate 
in space and grow exponentially in time. The electron plasma frequency for the number 
density Ue is given by cOp — {e^nej e^m^^l'^ and ^'(C) = — 2[1 + C-^(C)] is the first 
derivative of the plasma dispersion function = 7r~^/^ dt exp{—t'^)/{t — (). 

The plasma system is unstable for wave numbers in the range < k < kmax 
as depicted by the Fig. 1, where kmax — A^^^ depends on the temperature ratio 
T\\/T±_ > 1. The instability generates a magnetic field in perpendicular direction 
B = B± = (By + BlY^"^ and an electric field in parallel direction E = E\\ = E^. 
The connection between E and B is given by the Maxwell relation k x E = a;cB. 

2.2. The simulation method and parameters 

The PIC simulation is a method to model self-consistently the interplay of the electric 
and magnetic fields with the particles of a collision-less kinetic plasma [39] . The plasma 




(1) 




(2) 
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Figure 1. The growth rate (7(k). After a steep rise at low wave numbers the maximum 
is reached. Thermal effects limit the growth rate at large k. The wave number cutoff 
is given by A^/^. 



is treated as an incompressible phase space fluid, which is approximated by an ensemble 
of volume elements. The Vlasov equation is solved by the method of characteristics, 
i. e. a Lorentz equation of motion is solved for each volume element. A volume element 
contains a large number of physical particles, which are treated as one computational 
particle (CP) with the same charge to mass ratio as the physical particles. The velocity 
of a CP is at the position x^p. With the relativistic momentum p^p = f^cpl^cp 
that is given in terms of m^c and iricp in terms of the electron mass me, the normalized 
Maxwell equations for the electric and magnetic fields E(x, t) and B(x, t) 

VxE.-f.VxB.f + J, (3) 
and the normalized Lorentz equation 

X B[x,J) (4) 



dPcp 



dt 

are solved. The quantities E and B can be transformed to physical units by a 
multiplication by uOpiTLe/ce and uOpiTLe/e^ respectively. The normalized current density J 
is transformed to physical units by a multiplication with cngC. The qcp is given in units 
of the elementary charge e and the simulation time t and position x are given in units 
of l/uip and c/cUp, respectively. The code fulfills V • E = p/eo and V • B = to round-off 
precision [40]. 

The current J, the electric and the magnetic fields are defined on a grid. The 
electromagnetic fields are interpolated from the grid to the position of each CP. The 
momentum of the CP is updated with these interpolated fields and with the Lorentz 
equation of motion. The new particle position is calculated with the old particle position, 
Vcp and the time step A^. The microcurrents oc qcp^cp are the contributions of each 
discrete particle, i. e. the CPs in our simulations or the electrons in a real physical 
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plasma. They are interpolated from each CP back onto the grid and by integrating 
over all CPs the total current J is obtained. Then the electric and magnetic fields are 
updated with (3) and the individual steps are repeated. 

Our PIC simulations represent either one or two spatial dimensions and all 
components of Pcp- The boundary conditions are periodic in all directions. The 
quadratic 2D simulation box employs an equal number of cells — Ny — 1000 in 
the X- and in the y-direction, each with a side length A^ujp/c — 8.4 x 10~^. The number 
of CPs per cell is N2 — 150. For the ID simulation we use Ny — 1000 cells in the 
y-direction, N^; — 1 and the same A^;. The number of CPs per cell is in this case 
Ni = 32768. The time step is given by AtUp ~ 3.9 ■ 10~^. The equations solved by the 
PIC code arc given in normalized units and they can be scaled to any value of the plasma 
frequency. We choose a reference value for the plasma frequency {oUp — 6.3 • 10^ s~^) to 
illustrate the field amplitudes in physical units. 

3. Two-dimensional simulation 

3.1. The electric and magnetic energy densities 

The box-averaged energy densities of the electric Ei{x,y,t) and magnetic Bi{x,y,t) 
field components are given by eEi{t) = (TV^. A^^)"^ E eof-EjO'A^., A;A^, i)]^/2 and eBi{t) = 

{NxNy)~^Y,[Bi{jA.x,kA.x,t)f'/2iJiQ, respectively. The box-averaged kinetic energy 

density is ex(^) = {NxNy)~^A.~^Y,i^cpC^{.l3 — l); where the summation is over the 

3 

Ny Ni or NxNy N2 particles for the ID or the 2D simulation, respectively. The only 
magnetic field component that grows in the considered geometry in response to the 
TAWI is the component, which is confirmed by the PIC simulations. 

Fig. 2 demonstrates that the energy density of increases exponentially first and 
saturates for large simulation times similar to [22]. The electric field component Ey is 
also amphfied. The energy density of Ey does not grow in the 2D box at twice the 
exponential rate, with which the energy density of B^ is growing. It does so only in the 
ID simulation (shown here for comparison. Fig. 2(b)). The growth of the component 
E^ is delayed with respect to Ey. Eventually, both components E^ and Ey reach the 
same energy density in Fig. 2(a). In reference [22] the electric energy density is just 
a horizontal line because the energy densities have been plotted on a linear scale. The 
energy densities of the other field components did not grow during the simulation time. 

3.2. The field evolution 

Fig. 3 displays the magnetic field component cBz at the simulation times tuOp = 59 and 
tUp = 164. The left panel shows the magnetic field component cBz at the time when Ey 
saturates first. The component E^ has not been amplified yet. Magnetic filaments have 
developed with a symmetry axis parallel to the x-direction. The right panel displays 
the magnetic field component when both electric field components have reached almost 
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Figure 2. lO-logarithmic plot of the ratios eBz/£K{t — 0) (upper solid), eEy/^K{t — 0) 
(dashed) and eEx/^K{t — 0) (lower solid) for the 2D (a) and ID simulation (b). The 
TAWI amplifies B^, which induces a component Ey. In 2D is amplified later than 
Ey and their energies eventually equalize. 




Figure 3. The magnetic field component cB^ at the simulation times ttOp = 59 (a) 
and tiOp — 164 (b). The colourbar denotes the value in V/m (colours online). After the 
onset of the amplification of Ex the filaments are not strictly one-dimensional anymore. 



the same amplification level. The magnetic field structures are not strictly parallel to 
the a:-axis anymore (e.g. at yUp/c = 6.5), but they remain qualitatively unchanged in 
the considered reduced geometry. The structures would merge if the two perpendicular 
directions would be resolved [22]. 

The mechanism that transfers the electric field energy density from the Ey to the 
Ere component is revealed by comparing the spatial field profiles before and after the 
growth of Ex in the Fig. 4. The i?y-component is practically uniform along x and the 
structures are parallel to those of cBz at tujp = 75 and both are likely to be correlated. 
The amplitude Ey performs a full oscillation in the spatial interval 2.6 < yujp/c < 3.6 
between two zero-crossings of cBz- No significant structures can be seen in E^ at this 
time. The cB^ component is qualitatively unchanged in the time interval 75 < tUp < 87 
also on this scale, but the oscillation amplitude along y decreases somewhat. The electric 
field does, however, undergo a drastic change. The Ey that is initially uniform along x 
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Figure 4. The relevant field components at two times in a box sub-interval in units of 
V/m (colours online): The cBz, and Ey components at tujp — 75 are shown in the 
panels a,b,c, respectively. The panels d,e,f display the cB^, E^ and Ey components at 
ttUp = 87, respectively. 



becomes oscillatory at the later time. One oscillation takes place on a distance O.Sc/ujp 
along X, which corresponds to the wavenumber kxc/up ~ 8, as can be seen best from 
the Ey component. The oscillations along x rotate the electric field polarization vector 
away from y at tujp = 87. 

The observed spatial periodicity along x is evidence for a secondary instability, 
possibly a sausage mode instability [43]. It is this instability that couples the energy 
from the Ey component into the E^ component (see Fig. 2). This is confirmed by 
the supplementary movie 1, which animates in time the 10- logarithmic spatial power 
spectrum \E{k,j:, ky)\ of E{x,y) = Ex{x,y) + iEy{x,y). The movie 1 demonstrates that 
the electric field power is concentrated at A;^, ~ until tujp ~ 80, which is equivalent to 
a spatially uniform electric field along x. A signal grows after tUp ^ 80 at k^c/uip ~ 8, 
which corresponds to the oscillation of Ey in the Fig. 4(f). This signal is spread over 
a wide band of ky. Initially this signal grows only in the quadrant k^, ky > 0, but it 
rapidly spreads out in the wavenumber space after tUp 100. In [22] this secondary 
instability was not observed as only the magnetic field has been analyzed. 

3.3. The electron phase space distribution 

The variation of Bz along y must be supported by a current component Jx that varies 
along y. In order to obtain information about the distribution of the electrons in terms 
of the normalised momentum px/meC, we slice out 6 cells (0.084 < xuOp/c < 0.134) in 
the x-direction and integrate the phase space density over this x interval to improve the 
particle statistics. We compare the electron phase space distribution f{y,Px) to Jx{y)- 
Fig. 5 displays f{y,Px) and Jx{y) for the two simulation times corresponding to Fig. 
3. At tujp = 59 the mean momentum of the bulk electrons changes piecewise linearly 
with y and it spans a range {px/fnel < Vth\\- This results in a piecewise change of the 
current Jx- The B^iy) is consequently varying quadratic over extended y intervals (not 
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shown). At tUp = 164 these phase space structures become more diffuse. The magnetic 
field amphtude in Fig. 3 decreases by the reduction of at late times and, consequently, 
the energy density of in the Fig. 2(a). 




Figure 5. Top: The lO-logarithmic momentum distribution f{y,Px) at the simulation 
times tu!p = 59 (left) and tojp = 164 (right) (colours online): The structure of the 
electron phase space distribution is initially an approximately piecewise linear function 
of y. The structure is diffuse at late times. Bottom: Jx(y)l<lcp at tuip = 59 (left) 
and tujp = 164 (right). The current initially varies piecewise linearly with y. It is 
smoothened at late times and its amplitude has decreased. 



4. One-dimensional simulation 



We examine now with such a ID simulation the time beyond tUp ~ 50, when the TAWI 
grows the Ey field. The one- dimensional geometry suppresses E^. 




Z,„ ,4, 6 
(Panel a) yta^ic 




(Panel o) y a) / o 



Figure 6. The development of the field components cBz, Ey and during the 
simulation. The colourbar denotes the value in V/m. (a): The magnetic field structure 
is quasi-static and oscillatory in space, (b): Ey has a double frequency compared with 
cBz. Fast waves appear for tcjp > 180. (c): The low frequency waves of E^ are 
generated by Ampere's law with a resulting phase shift of 90° towards Bz- 



As in the 2D simulation only the z-component of the magnetic field has developed. 
The magnetic field structures are quasi-static and oscillatory in space (Fig. 6(a)). The 
structure of Ey in Fig. 6(b) has the half wavelength of cBz like in Fig. 3(a). For 
90 < tujp < 180 the structure is time- independent. This stationarity is an artifact of 
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the ID geometry. It permits us to separate the effects of and Ey on the electrons, 
which simpUfics their interpretation. For tUp > 180 fast waves appear in Ey, which are 
superposed on the periodic structure. Their velocity f ~ c/20 is of the order of Vth\\- 

The Ex in Fig. 6(c) resembles cB^ for tUp < 120 and has an amplitude that 
is one order of magnitude below that of Ey. Eq. (3), which is in one dimension 
dyBz — dtEx + J^, imphes that a growing generates the oscillations of E^. The phase 
shift of 90° between cB^ and E^ shown in Fig. 6 is an evidence for this connection. The 
oscillations of Er^ correspond to transient waves and they damp out as evidenced also 
by Fig. 2(b). 

Movie 2 shows the evolution of the electron charge density and the field components 
cB-^ and Ey. Initially, the growth of the cB^ and Ey components reveals a spatial 
correlation and the electron charge density is modulated accordingly. Oscillations of 
the electron charge density and of the Ey component with a short wavelength develop 
at late times. These rapid spatial oscillations obscure the long spatial oscillations. The 
latter are, however, still present, as the Fig. 6 demonstrates for Ey. The oscillation of 
B^ with the initial wavenumber remains dominant. 

In Fig. 7 snapshots are taken and the connection with the phase space distribution 
is investigated. In the ID simulation the electrons can move freely only along the y- 
direction. The TAWI generates initially a zigzag distribution, which yields through its 
currents the growth of the magnetic field. The amplitude of the velocity oscillations 
upon saturation are comparable to Vfh^^, which is clearly recognizable at tujp = 96 (See 
left column of Fig. 7 and movie 3). The wave number of Ey is twice that of cBz at 
this time and the zero crossings of Ey coincide with the extrema and zero crossings 
of cBz. At this stage the magnetic field still grows exponentially. Forces due to Ey 
and VcpB^ are not yet strong enough to affect the electron trajectories. Until strong 
electromagnetic fields have developed, the electrons perform just a thermal motion in 
the y-direction. Once the electromagnetic fields become strong at tcup — 116, they 
push together the particles along the ^/-direction, yielding a layered structure in the 
phase space distribution (middle column Fig. 7). The complexity of the phase space 
distribution increases further with the time, which is evident from the movie 3 and the 
right column in Fig. 7. The fine structure at the beginning of the simulation thermalizes 
but the zigzag distribution of the bulk remains. This holds up the currents which are 
neccessary to support the magnetic field. A coalescence of the filaments is not apparent. 

In order to demonstrate quantitatively the link between the magnetic pressure 
gradient and Ey, these are plotted in Fig. 8(a). The expected proportionahty of the 
magnetic pressure gradient and the electric field 

-eEy{y) = n:^ dyP,(y) = n:^i^^'B,(y) dyB,(y) (5) 

had to be modified by a factor 2. This factor can be understood as follows. The 
magnetic pressure gradient force accelerates the electrons, giving rise to a 7^ 0. The 
Jy is tied through Ampere's law to an Ey and the phases of both are shifted in time by 
90°. Initially, Ey — and Jy — and Ey must thus oscillate around an electric field 
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Figure 7. Row a: f{y,Px) of the particles at tUp = 96 (left), — 116 (middle) 
and = 154 (right). The colourbar denotes the lO-logarithm of the number of CPs. 

Row b: cB^ (dashed) and lOii'j, (solid) at the times corresponding to the electron 
distributions above. The y-axis denotes the value in V/m. 



set. The peak value of Ey is given by the added fields of the oscillatory Ey and of this 
set. This amounts to twice the amplitude of the set electric field, which is determined 
by the magnetic pressure gradient force. 




2 4,6 8 

y CO / c 



Figure 8. The normalized magnetic pressure gradient 2dyPb/nee (solid) and the 
negative electric field component —Ey (dashed) in values of V/m at tcop = 96. 



Fig. 9 analyses the electron distribution and the fields in a subsection of the 
simulation box at the time tcup — 96. The electron phase space distribution follows 
the zigzag distribution. The width of the distribution in with more than 300 particle 
counts is practically constant along y on the linear and 10-logarithmic density scales. 
The thermal spread and the electron density increase somewhat at the breaks, which 
coincide with the zero-crossings of Ey and ci?^. 

A comparison of the Figs. 9(a,b,d) reveals the cause of the density accumulation 
and the zigzag distribution. Let us consider the interval to the right of the break at 
yujp/c 2.1. There, the electron bulk has a p^; > and the Ey, < 0. The Ey < 
implies that the electrons accelerate away from the break in the direction of increasing 
y. The electrons accelerate away from the break also for the interval just to the left of 
yujp/c = 2.1, where Ey > 0. The break is thus an unstable equilibrium in what the Ey 
driven by VP5 is concerned and it cannot explain the electron accumulation close to 
yWp/c = 2.1 in Fig. 9(d). This is accomplished by the Bz- The Lorentz force component 
along y due to B^ for a computational electron with q^p < 0,Vcp,x > 0, B^ < just to 
the right oiyuip/c = 2.1 is F^p = —qcpVcp,xBz < and the electron is accelerated towards 
yujp/c — 2.1. The particles are also accelerated towards the break just to the left of it, 
because B^ changes its sign. It is thus the drift electric field oc v^B^ that compresses 
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Figure 9. The electron distribution and the electromagnetic fields in a sub-interval 
of the box at tuip — 96 (colours online), (a) The electron phase space density in units 
of a CP on a linear scale. The density decreases from the center outwards, (b) The 
electrostatic field lOEy (solid) and the cB^ (dashed) due to the TAWI. (c) The 10- 
logarithmic electron phase space density in units of a CP. (d) The electron charge 
density normalized to its mean value. 



the electrons. The electrostatic force that is due to the magnetic pressure gradient, 
which oscillates twice as fast as along y, opposes this compression close to the break 
and enforces it further away. The force contribution by the electrostatic field is not 
negligible. The peak speed along x of the electrons is 0.1c and the Lorentz force v x B 
is thus comparable to that due to Ey. 

5. Discussion 

In this work we have examined the thermal anisotropy-driven Weibel instability 
(TAWI). We have limited our investigation to an unmagnetized, non-relativistic plasma, 
considering only the electron dynamics and neglecting the ions. The TAWI is driven 
by a temperature anisotropy in the electron distribution. Here, the direction of higher 
temperature has been aligned with the x-direction. Thus, the TAWI modes would have 
wave vectors in the plane orthogonal to x. 

We have performed two simulations. The first one is a 2D simulation in the x-y- 
plane. As the electrons can move freely only in this plane, the waves that grow due 
to the TAWI are planar with a wave vector along y. Second we have performed a ID 
simulation along y with a large number of computational particles per cell. This has 
enabled through an increased dynamical range the identification of fine structure in 
the electron phase space distribution and it has reduced the noise. Both simulations 
have not resolved the y-^-plane, which would be necessary to enable the interplay of 
the current filaments [22]. We can thus separate the dynamics of individual current 
filaments driven by the TAWI from their mutual nonlinear interactions, as in [22], and 
leave the latter to future work. 

As expected, only the z-component of the magnetic energy density has been 
amplified [22]. After an initial exponential growth phase of Bz, the wave has saturated. 



PIC Simulations of the Weibel Instability 



13 



The growth of Bz has been accomphshed by the turning of the initially spatially uniform 
electron velocity distribution into a zigzag distribution in the phase space. The mean 
speed along x varies for this distribution piecewise linearly as a function of y. The initial 
Jx{y) = is thus transformed into a Jx{y) 7^ in most parts of the simulation box. In 
both simulations the saturation resulted in the growth of an electrostatic Ey component. 
We have demonstrated with the ID simulation the hnk between this Ey and the force due 
to the magnetic pressure gradient VP^. This connection was qualitatively unchanged 
in the 2D simulation, but the box-averaged electrostatic energy density did not grow 
at twice the exponential rate of the B^, as in the ID simulation of the TAWI and in 
simulations of the filamentation instability [18,41,44]. The E^ component could grow 
to a high amplitude only in the 2D simulation. The E^ in the ID simulation has been 
driven by the growing according to Ampere's law but only to a low amplitude. These 
modes damped out rapidly. The stronger E^ in the 2D simulation developed in response 
to a secondary instability, probably a sausage instability [43]. The ID simulation cannot 
resolve this instabihty, because it involves wavevectors along x. 

The ID PIC simulation provided an insight into the hnk between the fields and the 
electron phase space structures and it could resolve the charge density oscillations that 
developed during the nonhnear stage of the TAWI. The breaks in the zigzag distribution 
have been examined in more detail. The breaks correspond to the positions, at which the 
change with y of the mean momentum along x reverses its sign. The charge density and 
the electron temperature are increased at the breaks and the amplitudes of Ey, By have 
zero-crossings. The Ey field accelerates electrons away from the break, while the drift 
electric field imposed by the non-zero mean speed along x and the B^ is compressing 
the electrons along y. This compression is stronger close to the breaks than the electric 
repulsion, which results in an electron accumulation. 

We have demonstrated that the force due to the electric field is comparable to the 
Lorentz force imposed by the magnetic field, because the electron speeds are only a few 
per cent of the light speed. The focus on the magnetic forces and the neglection of 
the electric ones in Ref. [22] may thus not always be justified for the Weibel instability. 
The high phase space resolution of PIC simulations, which is now possible, has revealed 
that the size of the filaments in the ID simulation increases not through the merger 
of filaments, which is inhibited by the ID geometry, but through a broadening of the 
filaments due to the overlap of phase space layers. The filament mergers suggested by 
Fig. 3 in Ref. [22] could be an artifact of the lower plasma resolution or from the display 
of phase space slices instead of the full distribution. This has to be investigated with a 
future simulation that employs the anisotropy A = 24 of Ref. [22] . 

The connection of the magnetic and electric fields driven by the TAWI thus 
resembles those reported for the filamentation instability, which is driven by 
counterpropagating electron beams. Stable plasmons developed in [44] after the 
filamentation instability has saturated, which were confined by the potential due to 
the drift electric field and the magnetic pressure gradient. Here, the electron phase 
space distribution in the ID simulation developed fine structures first evolving into a 
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thermalized distribution with a bulk zigzag distribution. 

We will investigate in future work the statistical properties of the TAWI in ID and 
in 2D, as it has been done for the filamentation instability [42,45]. Ultimately, a 3D 
simulation has to be performed as in Ref. [23], which considered the original Weibel 
instability in which the electrons are hot along two cartesian directions and cool along 
the third. These simulations may provide an insight into the scale size and field strength 
of the magnetic field structures that develop out of nonthermal electron distributions 
in astrophysical fiows and to what extent these structures could be responsible for the 
frequently observed synchrotron radiation. 
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